Er 2 Ti 2 07: Evidence of Quantum Order by Disorder in a Frustrated Antiferromagnet 
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Er2Ti207 has been suggested to be a realization of the frustrated (111) XY pyrochlore lattice 
antiferromagnet, for which theory predicts fluctuation-induced symmetry breaking in a highly 
degenerate ground state manifold. We present a theoretical analysis of the classical model compared 
to neutron scattering experiments on the real material, both below and above Tn =1.173(2) K. The 
model correctly predicts the ordered magnetic structure, suggesting that the real system has order 
stabilized by zero-point quantum fluctuations that can be modelled by classical spin wave the- 
ory. However, the model fails to describe the excitations of the system, which show unusual features. 
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An important aspect of condensed matter is the sep- 
aration of energy scales, such that the minimization of 
one set of interactions may result in the frustration of 
another. A paradigm is the frustrated antiferromagnet, 
in which the local magnetic couplings between ions are 
frustrated by the crystal symmetry that the ions adopt. 
However, a systematic study of the rare earth pyrochlore 
titanates R,2Ti207 has shown that local antiferromag- 
netic bond frustration is neither a necessary, nor a suffi- 
cient condition for magnetic frustrationiaSiiLatti. Rather, 
it arises from the interplay, in the context of the crystal 
symmetry, of the principal terms in the spin Hamiltonian. 
In the case of R2Ti20y, the main terms are single-ion 
anisotropy, exchange and dipolar coupling. Depending 
on the balance of these factors, one observes spin ice be- 
havior (R = Ho, Dy)^, spin liquid behavior (R = Tb)^, 
and dipole induced partial order (R = Gd)£. 

Such behavior is best classified in terms of the domi- 
nant (111) single-ion anisotropy that arises from the trig- 
onal crystal electric field (CEF) at the rare earth site. For 
example, whereas the Heisenberg antiferromagnet has a 
spin liquid ground stated, the (111) Ising (dipolar) ferro- 
magnet has a spin ice ground stat o 1 ^ 7 . There is thus a 
clear motivation to study models based on other simple 
anisotropics and their realization in the titanate series. In 
this Letter we study one such model - the (111) XY model 
antiferromagnet 8 - and its realization Er2Ti2O*^ii iiiii£. 

We consider the Hamiltonian: 
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where the classical spins, Si, populate a face centered 



cubic array of corner sharing tetrahedra: the pyrochlore 
lattice. The spins are confined to easy XY planes by a 
local di — (111) anisotropy, D < 0, and are coupled an- 
tiferromagnctically by exchange J < 0. This model was 
first studied in Refj 8 ", where a discrete, but macroscopi- 
cally degenerate, set of ground states was identified. At 
finite temperature thermal fluctuations were found to se- 
lect an ordered state by the mechanism that Villain called 
"order by disorder"— and a first order phase transition 
was observed in numerical simulations. The propagation 
vector of the ordered state was found to be k = 0, 0, 
(henceforth "fc = 0"), but the basis vectors of the mag- 
netic structure were not determined. We have recently 
discovered that the ground state degeneracy is more ex- 
tensive than suggested in Refi 8 -; that there exists a con- 
tinuous manifold of k = ground states and that there 
may even be disordered states with continuous internal 
degrees of freedom— 

The possible basis states of the k — manifold were 
identified by group theory methods. They transform 
as linear combinations of the basis vectors of four irre- 
ducible representations (IRs), labelled ^5^9— The XY 
anisotropy energy is minimized only by (a) linear combi- 
nations of the two basis vectors tpi,2 that transform as the 
second order IR T 5 , (b) the discrete set of symmetrically 
equivalent basis vectors ^3-5 that transform as the third 
order IR (see Fig. 1). Monte Carlo simulations con- 
firmed that the first order transition at Tn/J = 0.125 se- 
lects k = order. An analysis of the distribution of bond 
energies suggested that immediately below Tn there re- 
mains a continuous degeneracy in the k — manifold 12 . 
As T — > 0, the spins were found to settle gradually into 
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FIG. 1: Upper inset: tetrahedral basis projected down [001], 
from left to right: ipi, ip2, ^3- ± denotes how spins tilt out 
of the plane (they lie parallel to the opposite triangular face 
of the tetrahedron) . Graph: size dependence of the simulated 
specific heat for N = 16L 3 spins. Lower inset: experimental 
low temperature specific heat vs T 3 (powder sample). A turn 
up low at temperature, attributed to hyperfine effects-, is 
barely visible with this choice of scale. 

the magnetic structure defined by 1^2 belonging to the 
continuously degenerate IR T 5 . (Note that ^2 is the only 
non-coplanar structure among 0i_5.) Both the initial 
selection of k = and the final selection of -02 must be 
order by disorder processes as the ground state manifold 
is macroscopically degenerate. 

To understand this ground state selection, we ana- 
lyze the D/J — > 00 model (i.e. spins confined to local 
XY planes). We expect that for the preferred ground 
state a spin wave analysis should expose the presence 
of zero frequency modes over an extensive region of the 
Brillouin zonei^. This is indeed the case: we calculate 
the quadratic Hamiltonian for small displacements away 
from a given ground state, which we symmetrize and di- 
agonalize to find the normal mode spectrum^. Applying 
this procedure to the state ip2 gives eigenvalues 

X(q) = 4J(l±cos(<f-a/2)) 

\(q) = 4j(l±cos(V-(c-6)/2)), (2) 

where a, b and c, are the basis vectors of the primitive 
rhombohcdral unit cell. Hence, there are branches with 
X(q) = over planes in the Brillouin zone, for which 
q- a = and q ■ (e— b) = 0. The same procedure when ap- 
plied to other selected ground states yields a microscopic 
number of zero modes at specific points in the zone. This 
difference gives the mechanism for the order by disorder 
selection of ip2- In this approximation the amplitude of 
the soft modes diverges, giving a dominant contribution 
to the entropjii^. Evidence of the soft modes exists in 
the specific heat, as each contributes less than ifcs. As 
there are 0(L 2 ) modes, the quantity \ — Ch/Nks should 
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FIG. 2: Main picture: powder neutron profile refinement 
(POLARIS, 50 mK, R wp = 1.23 %, about a third of the inten- 
sity is magnetic). The lower line is observed minus calculated 
intensity. Inset: single crystal Bragg intensity (PRISMA, 
(2,2,0) reflection) and fit to a power law (see text). 



scale as 1/L at low temperature, as confirmed in Fig. 1. 
We note that the entropy contribution to the free energy 
from the soft modes scales as Ns and so is not extensive. 
While this could mean that the ordering within the k = 
manifold occurs at a temperature-dependent system size 
that goes to zero in the thermodynamic limit, no such 
effect was detected in the system sizes we have studied. 

Disordered states are occasionally formed in the sim- 
ulations, by the rotation of columns of spins with infi- 
nite length out of an ordered state. This suggests that, 
starting from the ip2 state, one can introduce 0(L 2 ) in- 
dependent column defects, all perpendicular to a given 
plane. Our calculation, giving 0(L 2 ) soft modes is com- 
patible with this description and is analogous to the case 
of the Hciscnbcrg kagomc antiferromagnet 15 . For the lat- 
ter, fluctuations out of a coplanar spin configuration can 
be described equivalently in terms of soft propagating 
modes and localized zero energy excitations. 

The material Er2Ti2C>7, which orders magnetically at 
~ 1.2 KS, has been suggested to approximate the (111) 
XY antifcrromagnelpiSiii. To test this we have deter- 
mined its magnetic structure by powder neutron diffrac- 
tion using the POLARIS diffractometer (ISIS). The mag- 
netic reflections observed below Tn ~ 1.2 K index with 
a propagation vector of k = 0. As the transition is 
continuous (see below), the system is expected to order 
under only one of the non-zero IRs of the Er site rep- 
resentations: r3.5,7.9, as defined above. Refinement of 
the magnetic structure^ showed that only the two basis 
vectors tpi, ip2 of were consistent with the magnetic 
intensity (see Fig. 1, Fig. 2). Single crystal diffraction 
data collected on the instruments E2 (HMI, Berlin) and 
PRISMA (ISIS), allowed us to distinguish between the 
two structures. The measurements were performed on a 
(~ 8 mm 3 ) crystal at temperatures down to 0.13 K. In 
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order to suppress the formation of multidomains due to 
the cubic symmetry, a magnetic field was applied along 
the [1 1 0] direction. Below Tn, we found that a field 
of 0.5 T caused the (2,2,0) magnetic Bragg peak to in- 
crease from 260 counts to 500 counts, while the other 
peaks remained approximately unchanged. This increase 
in the (2,2,0) by a factor of 1.9 ± 0.2 is consistent with 
the formation of a monodomain of the ip2 ground state. 
We can conclude that the zero field ordering pattern is 
also described by ip 2 , in agreement with the theory. 

The ordered moment at 50 mK (where ordering is 
essentially complete) is 3.01 ± 0.05 [is per atom. A 
CEF analysis — , following RefiS 3 -, predicts a single-ion 
Kramers doublet ground state with the wavefunction 
-0.5428| - - 0.2384| - |) + 0.562S||) + 0.3876||) - 
0.426 1 ^} ■ This corresponds to moments of 3.8 fiB and 
0.12 hb perpendicular and parallel to (111). In an an- 
tiferromagnet, zero-point quantum fluctuations signifi- 
cantly reduce the ordered moment from the single-ion 
value. The observed moment of 3.01 /i_b, thus, agrees 
well with quasi-classical ordering of local XY-like mo- 
ments with a magnitude fixed by the CEF. 

The dominant perturbation to the model Hamilto- 
nian in Er 2 Ti 2 07 are dipole-dipole interactions™ and one 
might speculate that these are the true cause for the ip2 
ordering. In the Heisenberg pyrochlore antiferromagnet, 
provided the ratio of dipolar energy to near neighbor ex- 
change is less than a critical value, Jdd/Jnn < 5.7, the 
ground state is an energetically selected k = state: the 
tp3 basis (Fig. l) 2 ?i 21 . As this state is also a ground state 
for the (111) XY model, it is clear that in the present case 
it will be stabilized by dipolar interactions. We expect 
Er2Ti2C>7 to have Jdd/Jnn << 5.7, and thus -03 order, 
contrary to observation. With a moment of 3 hb the near 
neighbor dipolar interaction has a magnitude —0.32 K 
per spin for ip^ and +0.06 K per spin for ipi,2- It is there- 
fore a substantial fraction of Tn- However, zero point 
quantum fluctuations also stabilize the ground state of 
an antiferromagnet and will favor the softer ip2 stated. 

We are therefore drawn to the conclusion that the most 
likely cause of the observed ip2 order in Ei~2Ti2 07 is its 
stabilization by zero-point quantum fluctuations, the ef- 
fect of which is captured by our classical spin wave calcu- 
lation. However, the close agreement of experiment and 
theory does not extend to the excitations of the system, 
as we now discuss. 

In our model, the density of classical spin wave states, 
g(X) is a constant, a result that is incompatible with the 
experimental specific heat, which shows a T 3 dependence 
up to ~ 1 K (inset, Fig. 1). Furthermore, in contradic- 
tion with theory^, the ordering transition was found to be 
continuous within experimental precision. The evolution 
of the magnetic Bragg peaks near the critical tempera- 
ture Tn = 1.173(2) K obeys a standard power law with 
the critical exponent (3 = 0.33(2), characteristic of the 
3D-XY model (inset, Fig. 2)2^. This observation of a 
conventional (3 provides a counter example to the idea-^ 
that pyrochlore antiferromagnets represent a new univer- 
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FIG. 3: a) Powder inelastic spectrum (PRISMA, 1.8 K, |Q| 
= 1.1 A" 1 ). Inset: wave vector dependence of the excitation 
energies fitted to a cosine variation, b) Single crystal inelastic 
scattering (IN14, lOOmK, the strong elastic signal is Bragg 
scattering). Inset: Powder inelastic scattering (IRIS, 50 mK), 
integrated over 15 spectra from |Q| = 0.4 to 1 A -1 . 



sality class, analogous to Kawamura's chiral universality 
class for triangular lattice antiferromagnets^. 

To investigate the excitations further, the dynamics 
of Er2Ti2C>7 have been investigated by inelastic neutron 
scattering. Measurements were carried out on the spec- 
trometers PRISMA (ISIS), RITA (Ris0 National Labo- 
ratory) and IN14 (ILL). For a powder sample, excitongS^ 
were found at 6.3 and 7.3 meV (Fig. 3a). Their dispersion 
can be fitted to E(Q) = ^A 2 - 2a 2 J(Q)A with A = 
6.38(1) meV, a^/J{0) = 0.35(1) meV 1 / 2 for the lower 
level and A = 7.39(1) meV, a^/JiOj = 0.21(1) meV 1 / 2 
for the upper level. Here, J(Q) is the Fourier transformed 
exchange and a is related to the matrix elements that 
connect the ground and excited state a 26 i 27 . The values 
of A compare with the CEF predictions of 6.4 meV and 
8.8 meV— , while J(0)/k B is implied to be at least 1.5 K. 
Dispersed excitons are also observed in Tb2Ti207^. 

In the lower energy scale, data from an image-furnace 
grown single crystal (~ 250 mm 3 ) reveal a shoulder of 
scattering on the elastic line extending out to ~ 1.5 meV 
(Fig. 3b). This shoulder, which does not appear to obey a 
strong wave vector dependence (Fig. 3b) , gradually weak- 
ens upon heating above Tn- High resolution scans on 
PRISMA (T = 0.08 K) and on a powder sample on the 



4 



IRIS spectrometer at ISIS (inset, Fig. 3b) resolved the 
shoulder into a peak at 0.4 meV, with a width outside 
the instrumental resolution. 

One possible explanation, that this peak is a very low 
lying CEF level, as invoked by Blote et aim to explain the 
anomalously large Curie- Weiss temperature, 9 = —22 K, 
is not consistent with the predicted CEF scheme— A 
more natural explanation is that it is a weakly dispersed 
optical magnon mode. Evidence for a gapped magnon 
spectrum is provided by the IRIS data (inset, Fig. 3b), 
which is at the level of background between 0.03 meV 
(~ 0.3 K) and 0.2 meV (~ 2 K), indicating an absence of 
magnons in this range. However, there is a problem with 
this interpretation: the T 3 specific heat (inset, Fig. 1) 
requires a density of ungapped excited states that in- 
creases quadratically with energy up to at least ~ 1 K. 
One possible way out of this conundrum is to invoke the 
existence of a hidden branch of excitations that do not 
couple directly to neutrons. It is interesting to note that 
frustrated, quasi-two dimensional SrCr^Ga^-^Oig and 
(H 3 0)Fe3(S04) 2 (OII) 6 have T 2 specific heats despite the 
absence of conventional magnons 2 ^^. 

A final point of interest is that the other erbium py- 



rochlores Er 2 GaSb07 and Er2Sn207 apparently do not 
order down to 50 mK-^. It is possible that in this fam- 
ily there is a kind of quantum critical point where as 
a function of the interionic couplings, quasi-classical or- 
der gives way to the quantum fluctuations of the sin- 
gle ion doublets. Similar effects might also be rele- 
vant for Tb 2 Ti 2 0'A22, Yb 2 Ti 2 0722 and spinels such as 
ZnCrzO^M. 

In conclusion, the simplest theoretical model of 
Er 2 Ti 2 07, in which the single-ion moment is fixed by the 
CEF, the dipolar interaction is ignored and the exchange 
is treated quasi-classically, yields the correct ground state 
properties and furnishes strong evidence of quantum or- 
der by disorder. However, it fails to describe the un- 
usual excitation spectrum of the system and does not 
encompass the behavior of Er 2 GaSb 2 07 and Er 2 Sn 2 07. 
A detailed understanding of Er 2 Ti 2 07 thus leaves an in- 
triguing challenge for future research. 
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